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Abstract. The 68 000 J-band light curves of variable stars detected by the OGLE survey in the Large and Small Magellanic 
Clouds (MCs) are fitted by Fourier series, and also correlated with the DENIS and 2MASS all-sky release databases and with lists 
of spectroscopically confirmed M-, S- and C-stars. Lightcurves and the results of the lightcurve fitting (periods and amplitudes) 
and DENIS and 2MASS magnitudes are presented for 2277 M-,S-,C-stars in the MCs. The following aspects are discussed: the 
if-band period-luminosity relations for the spectroscopically confirmed AGB stars, period changes over a timespan of about 
17 years in a subset of about 400 LPVs, and candidate obscured AGB stars. The use of a sample of spectroscopically confirmed 
variables allows me to show specifically that almost all carbon stars are brighter than the tip of the RGB, and occupy sequences 
A+,B+,C and D. It is shown (for the LMC where there is a sufficient number of spectroscopically identified M-stars) that for 
sequences A+,B+,C the M-stars are on average fainter than the C-stars, as expected from an evolutionary point of view and 
previously observed in MC clusters. However, this is not so for sequence "D", suggesting that the origin of the so-called Long 
Secondary Periods is not related to an evolutionary effect. The fraction of objects that has a period on sequence "D" is also 
independent of chemical type. Three stars are identified that have been classified as oxygen-rich in the 1970s and carbon-rich 
in 1990s. Possibly they underwent a thermal pulse in the last 20 years, and dredged-up enough carbon to switch spectral type. 
The observations over almost two decades seem to suggest that up to 10% of AGB variables changed pulsation mode over that 
time span. More robust estimates will come from the ongoing and future (microlensing) photometric surveys. A sample of 570 
variable red objects ((J — K) > 2.0 or (J — K) > 4.0) is presented in which most stars are expected to be dust obscured AGB 
stars. Estimates are presented for cut-offs in (J — K) which should be applied to minimise dust obscuration in K, and based 
on this, C- and O-star AT-band PL-relations for large amplitude variables in the SMC and LMC are presented. 

Key words. Stars: AGB and post-AGB - Stars: carbon - Stars: variables: general - Magellanic Clouds 
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m ' 1 . Introduction confirmed and expanded upon by Noda et al. (2002), Lebzelter 

k> ! et al. (2002), Cioni et al. (2003), Ita et al. (2004a,b) and Kiss 

rN . In the course of the micro lensing surveys in the 1990's, the & Be dding (2003, 2004). These works differ in the source of 

Ctf monitoring of the Small and Large Magellanic Clouds has re- me varia bility data (MACHO, OGLE, EROS, MOA), area (SMC 

vealed an amazing number and variety of variable stars. A big or LMCX associated infrared data (Siding Spring 2.3m, DE- 

impact was felt and is being felt in many areas of variable NIS> 2masSi sirius)> and se i ec tion on pulsation amplitude or 

star research, like Cepheids and RR Lyrae stars. Also in the infrared colours 

area of Long Period Variables (LPVs) and AGB stars there ™ ., , . „ , T , 

° The present paper considers the OGLE data for both LMC 

has been remarkable progress. Wood et al. (1999) and Wood , T . , , ,.„„. . , ., _ T ^ , . • 

r b v ' and SMC. Ita et al. (2004a) only consider the OGLE data in 

(2000) were the first to identify and label different sequences , ... . TT , . . T , ,„ , „, ,„ 

„._,„„ , , overlap with their SIRIUS IR observations in LMC and SMC, 

ABC thought to represent the classical Mira sequence ( C ) , „. „ ... ., . . 

° , .... „,„ . „^„ .... .,, and Kiss & Bedding consider only stars in the SMC with 

and overtone pulsators ( A,B ), and sequence D which is still _ , ( T K\ no 

unexplained (Olivier & Wood 2003, Wood 2003, 2004). Stars S } ^ 

on these sequence are sometimes referred to as having Long Also in contrast t0 P revious studies ' em P hasis is P ut on 

Secondary Periods-LSPs. This view has subsequently been spectroscopically confirmed AGB stars (i.e. M-, S- and C- 

stars). In other studies M- and C-stars are usually identified 

photometrically by using a division at a (J — K) colour of 

Send offprint requests to: Martin Groenewegen . . ™ . ... .. , , il „ 

_ , , , . ^1.4 mag. This paper specincally addresses the properties of 

(groen@ster.kuleuven.ac.be) ... u»t.,4™. 

i ui ■ i c . .u known carbon stars in relation to sequences ABCD . 

Tables 2, 3, 6, 7 and 8 are available in electronic form at the n 

CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via The paper is structured as follows. In Section 2 the OGLE, 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-i-A/ Figures 2, 9 and 13 2MASS and DENIS surveys are described. In Section 3 the 

are available in the on-line edition of A&A. model is presented, both in terms of the actual lightcurve fit- 
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ting, and the post-processing. Section 4 presents the results. 
Discussed are the /C-band PL-relation for the spectroscopi- 
cally confirmed AGB stars, period changes over a timespan 
of almost 2 decades, and a sample of very red obscured AGB 
star candidates. The conclusions are summarised in Section 5. 
Some of this work, and the star-to-star comparison of periods 
derived by me from MACHO and OGLE data and literature val- 
ues are described in Groenewegen (2004). 



2. The data sets 



Zebruri et al. (2001) describe the dataset of the about 68 000 
variable objects detected by OGLE in the direction of the LMC 
and SMC, obtained in the course of the OGLE-II micro lensing 
survey using the difference image analysis (DIA) technique. 
Twenty-one fields in the central parts of the LMC, and 1 1 fields 
in the central parts of the SMC of size 14.2' x57' each were 
observed in BVI, with an absolute photometric accuracy of 
0.01-0.02 mag. The large majority of data was obtained in the 
/-band (and the DIA analysis has been done only on the I- 
band data), and these data were downloaded from the OGLE 
homepage (http://sirius.astrouw.edu.pl/~ogle/ 1. 

The DENIS survey is a survey of the southern hemisphere 
in IJK S (Epchtein et al. 1999). Cioni et al. (2000) describe 
a point source catalog of sources in the direction of the 
Magellanic Clouds (MCs; DCMC = DENIS Catalogue towards 
the Magellanic Clouds) containing 1 635 680 objects with a de- 
tection in at least 2 of the 3 photometric bands that is available 
in electronic form. The 68193 OGLE objects were correlated 
on position using a 3" search radius and 40793 matches were 
found. 

The 2MASS survey is an all-sky survey in the JHK S near- 
infrared bands. On March 25, 2003 the 2MASS team released 
the all-sky point source catalog. The easiest way to check if a 
star is included in the 2MASS database is by uplinking a source 
table with coordinates to the 2MASS homepage. Such a table 
was prepared for the 68 193 OGLE objects and correlated on po- 
sition using a 3" radius. Data on 50129 objects were returned. 

3. The model 

At the heart of the data processing are two numerical codes, 
that are described in detail in the Appendices. Briefly, the first 
code (see for details Appendix A) sequentially reads in the I- 
band data for the 68 000 objects, determines periods through 
Fourier analysis, fits sine and cosine functions to the light curve 
through linear least-squares fitting and makes the final correla- 
tion with the pre-prepared DENIS and 2MASS source lists. All 
the relevant output quantities are written to file. 

This file is read in by the second code (see for details 
Appendix B). A further selection may be applied (typically on 
period, amplitude and mean /-magnitude), multiple entries are 
filtered out (i.e. objects that appear in different OGLE fields), 
and a correlation is made with pre-prepared lists of known non- 
LPVs and known LPVs or AGB stars. The output of the second 
code is a list with LPV candidates. 

The third step (for details see Appendix C) consists of a 
visual inspection of the fits to the light curves of the candi- 



figl.jpg 



Fig. 1. Difference in photometry, after the following offsets 
have been applied: /(denis-ogle) = —0.018, J(denis-2mass)= 
—0.090, A"(denis-2mass)= —0.14. Used in the analysis are the 
stars with an /-band amplitude < 0.05 mag. Plotted are about 
12000, 8700, and 5700 stars, top to bottom. 



date LPVs and a literature study through a correlation with the 
SIMBAD database. Non-LPVs are removed, and sometimes the 
fitting is redone. The final list of LPV candidates is compiled. 



4. Results 

4. 1. Astrometry 

The spatial correlation between the OGLE objects and known 
LPVs and AGB stars, and known non-LPVs, is actually done 
in 2 steps. In the first step the correlation is made, and the dif- 
ferences and spread in ARA cos(5) and A5 are determined. 
These mean offsets are then applied to make the final cross- 
correlation. The results are listed in Table ^ Regarding the 
MACHO data the agreement between the offsets determined 
from the list of non-LPVs, and the data by Wood et al. (1999) 
is good and a combined offset of A RA cos((5) = 0.52" and 
AS = 0.16" is applied. With respect to the other sources of as- 
trometry, similar small changes have been found, and in most 
cases applied, usually increasing the number of matches in the 
final matching. 
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Table 1. Comparison of coordinates, and number of positional matches, before and after a correction was applied. 



OGLE- (other) 


A RA cos (<5) 


AS 


N 


A RA cos (<5) 


A<5 


N 


remark 


MACHO 


0.52 ± 


0.61 


0.21 ± 0.63 


647 


0.06 ± 


0.63 


0.05 ± 0.64 


662 


mainly "blue objects" and EBs 


MACHO 


0.52 ± 


1.28 


0.14 ± 1.23 


1427 


0.11 ± 


1.31 


-0.02± 1.22 


1436 


LPVs; Wood et al. 1999 


MOA 


— 0.01± 


0.32 


0.008 ± 0.55 


123 










EBs; no correction applied 


MOA 


n 7x + 

U. / O I 


1 04. 


0.002 ± 0.94 


150 


0.11 ± 


n 07 

\j.y i 


0.03 ± 0.94 


150 


LPVs; Noda et al. 2002 


EROS 


0.41 ± 


0.50 


—0 21+ 48 


368 


0.06 ± 


0.58 


01 + 49 


381 




AGAPEROS 


0.31 ± 


1.09 


0.37 ± 1.00 


503 


0.05 ± 


1.07 


0.05 ± 0.95 


501 


LPVs; Letzelter et al. 2002 


DENIS 


-0.10± 


1.33 


0.02 ± 1.37 


313 










LPVs; Cioni et al. 2003; no correction applied 


HUGHES 


2.19 ± 


1.50 


0.85 ± 1.41 


433 


0.22 ± 


1.34 


-0.01± 1.34 


436 


LPVs 


KDM 


0.55 ± 


1.12 


-0.23± 0.90 


1024 


0.06 ± 


1.11 


-0.03± 0.90 


1028 


carbon stars 


RAW 


-1.06± 


0.80 


1.03 ± 0.71 


873 


0.04 ± 


0.76 


0.04 ± 0.70 


874 


carbon stars 


BMB 


0.42 ± 


0.90 


0.00 ± 0.66 


381 


0.05 ± 


0.96 


-0.01± 0.68 


385 


AGB stars 



Table 2. First entries in the electronically available table, which lists: OGLE-field, OGLE-name, the three fitted periods with 
errors and amplitude (0.00 means no fit), mean logic, and associated DENIS UK photometry with errors, and associated 2MASS 
JHK photometry with errors (99.9 and 9.99 means no association, or no value). 



Field 



Period 

95.2 ± 0.138 

57.8 ± 0.054 
0.0 ± 0.000 

83.9 ± 0.200 
135.3 ±0.1 18 
113.8 ±0.083 

179.2 ± 0.057 
107.5 ± 0.105 
188.1 ±0.377 

616.8 ± 12.389 
105.1 ± 0.040 

66.3 ± 0.070 
251.1 ± 0.970 
118.1 ±0.064 

0.0 ± 0.000 
95.3 ± 0.049 
361.0 ± 1.652 

172.3 ±0.218 
65.2 ±0.102 

102.0 ±0.101 
646.7 ± 6.002 

120.1 ± 0.150 

171.2 ± 0.486 



Amp 



Period 
168.5 ±0.335 
88.5 ± 0. 142 

0.0 ± 0.000 
166.9 ± 0.393 
755.8 ± 1.033 
186.5 ± 0.420 

0.0 ± 0.000 
135.3 ±0.198 
281.8 ±0.657 

0.0 ± 0.000 

251.8 ±0.291 
0.0 ± 0.000 
0.0 ± 0.000 

222.2 ± 0.1 IS 
0.0 ± 0.000 

189.9 ± 0.297 

481.1 ± 1.764 
216.8 ±0.170 
603.8 ± 3.661 

342.2 ± 0.527 
0.0 ± 0.000 

131.3 ±0.158 
243.0 ± 0.773 



Amp 



771.7 ± 7.449 

730.4 ± 4.627 
0.0 ± 0.000 
0.0 ± 0.000 
0.0 ± 0.000 

727.8 ± 3.174 
0.0 ± 0.000 

1330.7 ± 14.439 
0.0 ± 0.000 
0.0 ± 0.000 
559.2 ± 1.585 
0.0 ± 0.000 
0.0 ± 0.000 
0.0 ± 0.000 
0.0 ± 0.000 

514.9 ±2.064 

872.6 ±2.103 
0.0 ± 0.000 
0.0 ± 0.000 

684.9 ± 0.902 
0.0 ± 0.000 

328.7 ± 0.541 

542.5 ±6.147 



Period Amp I a , 



a;. 



■J'2\ I ASS 



II; 



'j\[.\SS 



K 



2.\[.\SS 



smcjscl 
smc_scl 
smc_scl 
smcjscl 
smcjscl 
smc_scl 
smc_scl 
smcjscl 
smcjscl 
smcjscl 
smc_scl 
smc_scl 
smcjscl 
smcjscl 
smc_scl 
smc_scl 
smcjcl 
smc_scl 
smc_scl 
sme_scl 
smcjcl 
smc_scl 
smc_scl 



OGLE003616 
OGLE003626 
OGLE003630 
OGLE003633 
OGLE003633 
OGLE003635 
OGLE003638 
OGLE003643 
OGLE003646 
OGLE003649 
OGLE003650 
OGLE003651 
OGLE003653 
OGLE003700 
OGLE003703 
OGLE003705 
OGLE003712 
OGLE003716 
OGLE003718 
OGLE003722 
OGLE003726 
OGLE003734 
OGLE003735 



76-732132.7 
98-733704.3 
80-735456.4 

16- 730547.6 
30-735040.6 
63-733551.9 
83-734635.4 
09-735018.3 
48-734220.6 
80-730416.4 
09-731737.3 
78-732400.0 

17- 733620.3 
95-731637.7 
93-735608.3 
38-734426.3 
83-731133.5 
26-732048.9 
97-735029.8 
15-732222.8 
67-731710.4 
50-732110.6 
61-730421.6 



0.042 
0.020 
0.000 
0.029 
0.063 
0.050 
0.243 
0.051 
0.035 
0.011 
0.090 
0.018 
0.021 
0.065 
0.000 
0.061 
0.104 
0.052 
0.011 
0.033 
0.023 
0.056 
0.023 



0.064 
0.018 
0.000 
0.057 
0.195 
0.027 
0.000 
0.043 
0.046 
0.000 
0.072 
0.000 
0.000 
0.126 
0.000 
0.043 
0.148 
0.098 
0.032 
0.100 
0.000 
0.063 
0.028 



0.140 
0.041 
0.000 
0.000 
0.000 
0.049 
0.000 
0.102 
0.000 
0.000 
0.070 
0.000 
0.000 
0.000 
0.000 
0.056 
0.315 
0.000 
0.000 
0.274 
0.000 
0.063 
0.019 



14.47 
14.73 
14.44 
14.24 
14.56 
14.22 
15.16 
14.21 
14.30 
14.52 
14.14 
14.53 
14.15 
14.42 
15.28 
14.72 
14.27 
14.39 
14.26 
14.36 
14.71 
14.61 
14.19 



14.42 ± 0.01 
14.86 ± 0.02 
14.54 ± 0.01 
14.23 ± 0.01 
14.47 ± 0.01 
14.31 ±0.01 
14.85 ± 0.02 
14.07 ± 0.01 
14.22 ± 0.01 
14.66 ± 0.01 
14.19 ±0.01 
14.62 ± 0.01 
14.22 ± 0.01 
14.51 ±0.01 
99.00 ± 9.99 
14.78 ± 0.01 
14.13 ± 0.01 
14.66 ± 0.01 
14.29 ± 0.01 
14.77 ± 0.01 
14.75 ± 0.01 
99.00 ± 9.99 
99.00 ± 9.99 



12.66 ± 0.01 
13.38 ±0.02 
13.09 ± 0.01 
12.68 ±0.01 
12.89 ± 0.01 
12.59 ± 0.01 

13.13 ±0.02 
12.43 ±0.01 
12.76 ±0.01 

13.14 ±0.01 
12.51 ±0.01 
13.04 ±0.01 
12.65 ± 0.01 
12.84 ±0.01 
99.00 ± 9.99 
13.22 ±0.01 
12.42 ±0.01 
13.06 ± 0.02 
12.98 ±0.01 
13.12 ±0.02 
13.24 ±0.02 
99.00 ± 9.99 
99.00 ± 9.99 



11.17 ±0.02 
12.32 ±0.06 

11.93 ±0.05 
11.48 ±0.03 
11.39 ±0.03 
11.06 ±0.02 
11.66 ±0.04 
11.12 ±0.02 
11.68 ±0.04 

11.94 ±0.05 
11.28 ±0.03 
11.78 ±0.04 
11.36 ±0.03 
11.45 ±0.03 
99.00 ± 9.99 
11.56 ±0.03 
10.81 ±0.02 
11.62 ±0.04 
11.81 ±0.04 
11.60 ±0.03 
12.15 ±0.06 
99.00 ± 9.99 
99.00 ± 9.99 



12.79 ± 0.03 
13.48 ± 0.03 
13.20 ± 0.03 
12.78 ± 0.02 
13.00 ± 0.03 
12.72 ± 0.03 
13.26 ± 0.03 
12.58 ± 0.02 
12.95 ± 0.03 
13.22 ± 0.03 
12.71 ± 0.03 
13.17 ±0.03 
12.75 ± 0.02 

13.06 ± 0.03 
14.20 ± 0.04 
13.28 ± 0.04 
12.53 ± 0.03 

13.09 ± 0.02 

13.07 ± 0.02 
13.28 ± 0.02 
13.35 ±0.03 

13.10 ±0.03 
12.71 ±0.03 



11.74 ±0.03 

12.72 ± 0.03 
12.38 ± 0.03 
11.90 ±0.03 
12.02 ± 0.03 

11.73 ±0.03 

12.28 ± 0.03 
11.67 ±0.03 
12.13 ±0.03 
12.41 ± 0.03 
11.80 ±0.03 

12.29 ± 0.03 
11.87 ±0.03 
12.06 ± 0.03 

13.44 ±0.05 
99.00 ± 9.99 

11.45 ±0.03 
12.15 ±0.02 
12.26 ± 0.02 
12.23 ± 0.02 
12.51 ±0.03 
12.08 ± 0.02 
11.73 ±0.02 



11.27 ±0.03 
12.34 ± 0.03 
12.12 ±0.03 
11.55 ±0.03 
11.52 ±0.03 
11.18 ±0.02 
11.74 ±0.03 
11.27 ±0.03 
11.81 ±0.03 
12.12 ±0.03 
11.46 ±0.02 
12.00 ± 0.03 
11.48 ±0.02 
11.59 ±0.03 

13.23 ±0.06 
99.00 ± 9.99 
10.94 ± 0.03 
11.76 ±0.03 
11.99 ±0.02 
11.76 ±0.03 

12.24 ± 0.02 
11.66 ±0.02 
11.32 ±0.02 



Table 3. First entries in the electronically available table, 
which list: OGLE-field, OGLE-name, other names, spectral 
type and references and comments. 



Field 


OGLE-name 


Other Name 


SpT Comments 


smc_sc 1 


OGLE003616.76-732132.7 


MH.0344 


C 


smc_sc 1 


OGLE003626.98-733704.3 


MH.0349 


C 


smc_sc 1 


OGLE003630.80-735456.4 


MH.0352 


c 


smc_sc 1 


OGLE003633. 16-730547.6 


RAW.0003 


c 


smcscl 


OGLE003633.30-735040.6 


MH.0354 


c 


smc_sc 1 


OGLE003635.63-733551.9 


RAW.0004 


c 


smc_sc 1 


OGLE003638.83-734635.4 


MH.0358 


c 


smc_sc 1 


OGLE003643.09-735018.3 


MH.0361 


c 


smc_sc 1 


OGLE003646.48-734220.6 


RAW.0005 


c 


smc_sc 1 


OGLE003649.80-730416.4 


RAW.0007 


c 


smc_sc 1 


OGLE003650.09-731737.3 


RAW.0008 


c 


smcscl 


OGLE00365 1.78-732400.0 


RAW.0009 


c 


smc_sc 1 


OGLE003653. 17-733620.3 


RAW.0010 


c 


smc_sc 1 


OGLE003700.95-731637.7 


RAW.0011 


c 


smc_sc 1 


OGLE003703.93-735608.3 


MH.0370 


c 


smc_sc 1 


OGLE003705.38-734426.3 


MH.0371 


c 


smc_sc 1 


OGLE003712.83-731133.5 


RAW.0012 


c 


smc_sc 1 


OGLE003716.26-732048.9 


RAW.0013 


c 


smcscl 


OGLE0037 1 8.97-735029.8 


MH.0374 


c 


smcscl 


OGLE003722. 15-732222.8 


RAW.0014 


c 


smcsc 1 


OGLE003726.67-731710.4 


RAW.0015 


c 


smcsc 1 


OGLE003734.50-732110.6 


RAW.0016 


c 


smcsc 1 


OGLE003735.61-730421.6 


RAW.0017 


c 



4.2. Photometry 

As was already discussed in Groenewegen (2004), by selecting 
the least variable stars in the OGLE database one can compare 
photometry. 

In particular the OGLE / was compared to the (singe-epoch) 
DENIS /, and the (single-epoch) DENIS JK was compared to 
the (single-epoch) 2MASS JK magnitudes. This was done by 
selecting those objects with an amplitude in the /-band of < 
0.05 mag. 

Figure[2shows the final results when offsets /(denis-ogle) 
= —0.018, J(denis-2mass)= —0.090, and i4'(denis-2mass)= 
—0.14 are applied. The offsets derived here are very similar 
to those derived in Delmotte et al. (2002) using a similar anal- 
ysis based on a direct comparison of DENIS with the 2MASS 
2nd incremental data-release (they found: J(denis-2mass) = 
-0.11 ± 0.06, and if (denis-2mass)= -0.14 ± 0.06). 

4.3. The sample of spectroscopically confirmed 
M,S,C-stars 

As mentioned before, a correlation was made with a list that 
contains 12631 objects over the entire MCs with a spectroscop- 
ically determined spectral type. Using a 4" search radius, 2478 
unique objects were found in the OGLE fields. After visual in- 
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spection (removing mostly obviously incorrect positional mis- 
matches, indicated by association with very faint OGLE objects 
with hardly any variability [sect. 3.2.2]) there remain 2277 ob- 
jects (856 C- and 3 M-stars in the SMC, 1064 C-, 10 S-, 344 
M-stars in the LMC) which are the subject of further study. 
Their lightcurves are shown in Figure [2] (only the first few, the 
full figure only being available in the electronic edition), the 
results of the lightcurve analysis and association with DENIS 
and 2MASS is given in Tabled and the association with known 
objects and additional comments and references are given in 
Table [5] (again, only the first few entries are shown, the full 
tables being available only electronically). 

In the discussion that follows, magnitudes are dered- 
dened using the Ay values that correspond to the respective 
OGLE field in the SMC or LMC, and selective reddenings of 
Aj/Ay = 0.49, Aj/Ay = 0.27, A n /Ay = 0.20, A K /Ay = 
0.12 (Draine 2003) are used. 

Some of the objects have different spectral classifications 
in different surveys. These have been identified in Table [3] 
Some of them are in fact due to a typographical error in CML 
(Cioni, private communication). In all other cases the differ- 
ent sources of photometry, pulsation periods (when available) 
and close proximity suggest that these objects are one and the 
same. For one object, OGLE05271 1.00-692827.5, both associ- 
ated objects have an uncertain spectral type: WBP-48 is clas- 
sified as (M?), SHV052733.4-693050 as (CS?). It is kept as 
"M". OGLE052705.07-693606.9 is associated with the car- 
bon star KDM-4094, but also with the M-star HV 12048, 
and the latter type is assumed here. Of course oxygen-rich 
stars may evolve into carbon-rich stars, and three more ob- 
jects that may have evolved in this way judging from the dif- 
ferent spectral types are OGLE050859.83-691458.0(WOH-G- 
202, KDM-2226), OGLE054109.00-700942.1 (WOH-G-473, 
KDM-5626), and OGLE054120.38-700823.3 (WOH-G-478, 
KDM-5645). In all cases the more recent spectral type (C) has 
been adopted. These stars are valid targets for spectroscopy, 
both in the optical to determine their (s-process) abundances, 
and in the mid-infrared to identify the dust features of the pos- 
sibly still oxygen-rich dust shell. 

Six objects are erroneously associated by SIMBAD with 
relatively bright objects, and an appropriate comment is added 
in Table [5] The alledged counterparts have B and/or V magni- 
tudes in the range 10-13 mag, while the / og i e , and the photo- 
electric R and / magnitude (all six objects happen to be listed 
in KDM) are typically 14-15 mag and suggest B, V magni- 
tudes that would be closer to 16-18 mag. Optical and 2MASS 
finding charts were also inspected and a bright object was typ- 
ically found within 1' of the OGLE object. 

One of the major finding by Wood et al. (1999), Wood 
(2000) and subsequently confirmed by other studies is the ex- 
istence of distinct PL-relations, usually indicated by the letters 
A-,A+,B-,B+,C,D. The values for the boundaries of these re- 
gions were originally taken from Ita et al. (2004a), Ita (private 
communication) and then slightly adapted. They are basically 
drawn by eye, and for reference the lines that define the bound- 
aries have been listed in Table |4] The boundaries between A- 
and A+, and B- and B+, are placed at the tip of the RGB, i.e. 
K = 12.1, and 12.7 in the LMC and SMC, respectively. 



Figure [3] shows the PL-relation in the LT-band with these 
boundaries, for the SMC and LMC separately, and according 
to cuts in the /-band amplitude, as indicated in the insets. The 
/T-band is on the 2MASS system, and is the average of the 
DENIS and 2MASS photometry. In particular, if both DENIS 
and 2MASS X-band data is available, the DENIS data point is 
corrected as explained above (i.e. 0.14 mag added), and aver- 
aged with the 2MASS data point. This should take out some 
of the scatter in the PL-diagram, as the effect of the vari- 
ability in the A'-band is reduced. If only DENIS is available, 
the corrected value is used. Not all periods listed in Table |2] 
are plotted. Inspecting the outliers in preliminary versions of 
Figure [3] indicated that to reduce the scatter a further culling 
was necessary, in particular at the longer periods, as the times- 
pan of the OGLE observations is about 1 100 days. The follow- 
ing conditions were applied (with AP the error in the period): 
AP/P < 0.01 for P < 500 d ; AP < b d for 500 d < P < 800 d 
and AP < 1.5 d for P > 800 d . Figures EJ E| 13 El EI CQl C21 ESI 
and all calculations are based on periods that obey these condi- 
tions. In contrast, all periods found have been listed in Tables|2] 
and0and are shown in Figures El|9]and[T3] 

The following systematics may be observed: (1) small am- 
plitude variables are present at all periods, (2) objects with I- 
band amplitudes ^0.05 are not found in box "A+" (nor "A- 
"), (3) objects with /-band amplitudes ,>0.45 are mainly found 
in box "C". These three remarks are valid for SMC and LMC 
alike. 

Unfortunately, only three confirmed oxygen-rich SMC 
stars appear in Tables [2] and [3] and all come from an IRAS 
selected sample (GB98). They all have (very) long periods and 
the error in these periods are such that only one period appears 
in Figure|3] 

For the LMC the situation is better. One can observe that: 
(4) essentially all SMC, and almost all LMC, carbon stars are 
brighter than the tip of the RGB, (5) LMC M-stars are observed 
below the TRGB, but only at small amplitudes, (6) for a given 
"box" or amplitude cut, the M-stars are on average fainter than 
the C-stars, (7) a few stars brighter than the expected tip of the 
AGB, presumably supergiants, are present, and they predomi- 
nantly have large pulsation amplitudes. 

Statement 6 is illustrated in more detail in Figures 0] and 
|5j where, respectively, LMC and SMC K -band histograms are 
shown for the stars inside different "boxes" and partly for dif- 
ferent /-band amplitude cuts, as indicated in the insets. For the 
LMC a distinction is made between C- and M/S-stars and for 
sequences A+,B+ and C the M-stars are on average fainter than 
the C-stars. The seems not to be the case for sequence D, sug- 
gesting that the LSP phenomenon is not related to evolutionary 
status. 

Figure [6] shows the Amplitude-Period diagram for LMC 
and SMC. They look qualitatively similar, and similar to di- 
agrams shown in Ita et al. (2004a) and Cioni et al. (2003). 
Interestingly, it appears that for a given amplitude, the LMC 
M-stars have on average shorter periods than the C-stars. This 
is a different manifestation of something that could already be 
seen in Figure|5J most clearly in the panel with the cut 0.15 < 
Amplitude < 0.45 where sequence "B+" contains mainly M- 
stars while sequence "C" contains mainly C-stars. 



M.A.T. Groenewegen: Long Period Variables in the Magellanic Clouds: OGLE + 2MASS + DENIS 



5 



fig2.jpg 



Fig. 2. First entries of electronically available figure with all lightcurves. The fit is indicated by the (red) solid line. Crosses 
indicate data points not included in the fit. 

Table 4. Definition of the boundary lines of the sequences in the K-bmd PL-diagram in Figure Relations are of the form 
-Kmin.max = slope m i n max log P + zp m i n max , with K the dereddened magnitude on the 2MASS system. 
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Fig. 3. K -band Pi-relation, for the SMC (left) and LMC (right). Panles indicate selection on /-band amplitude. Carbon stars 
are indicated by filled circles, M- and S-stars by open circels. All periods from Table|2]that fulfil AP/ P < 0.01 for P < 500 d ; 
AP < b d for 500 d < P < 800 d and AP < 1.5 d for P > 800 d are plotted. Boxes related to the "ABCD" sequences are 
indicated. The functional dependence is summarised in Table|4] 



Figure shows the Period-Colour diagram for LMC and 
SMC for DENIS and 2MASS data. They look qualitatively sim- 
ilar to diagrams shown in Lebzelter et al. (2002) and Ita et al. 
(2004a). Two features may be remarked upon, (1) there are 
more red carbon stars in the LMC than in the SMC (at least that 
are spectroscopically confirmed), (2) The M-stars are much 
more spread in (I — J)q than in (J — K)q colour. 

The first observation may be related to an, on average, 
higher mass loss rate of AGB stars in the LMC compared to 
the SMC. Section 1431 deals in more detail with obscured ob- 
jects. The second observation is a known effect and related to 



the strong temperature dependence of the molecular bands in 
the red part of the optical spectrum in M-stars. Table|5]rists the 
(I — J) and (J — K) colours according to the model atmo- 
spheres (for solar metallicity) of M-stars (Fluks et al. 1994), 
and carbon stars (Loidl et al. 2001, her model with a C/O ra- 
tio of 1.1). In addition, for a few selected models the colours 
are listed in the last column of the stars when obscured by 0.05 
magnitudes in K by dust (see Section l4~5Y 




Fig. 3. Continued 



4.4. Changes in periods over time 

In this section possible changes in pulsation periods over time 
are being discussed. This is possible because the periods in 
the Hughes list and WBP are based on plate material typically 
taken between 1977 and 1984, while the MACHO, OGLE, 
MOA, AGAPEROS data was taken in the late-1990s, i.e. a 
timespan of typically 17 years. A cross-correlation of the OGLE 
objects was performed, but this time not with spectroscopically 
known objects, but rather with known LPVs with a historical 
period. Table [6] (complete table available in electronic form) 
lists the results for a total of 370 stars and is similar to Table|2] 
In addition the known periods are listed: first the historical pe- 
riod (Hughes list or WBP), then the corresponding period from 
OGLE, and when available other known periods (with reference 



between parenthesis). Since multiple periods are allowed for by 
the analysis of the OGLE data, the OGLE period quoted is the 
one for which the amplitude is similar to the one corresponding 
to the historical period (in some cases this then corresponds not 
to the primary period as determined from the OGLE data, but 
often to the LSP). In the last column some remarks are given 
concerning possible period or pulsation mode changes. 

Of particular interest to look for period changes are the 
Miras that have defined the Mira Pi-relation over the years 
and that have periods determined at multiple epochs (GLE81, 
WBP, FGWC, GLE03). In fact, GLE03 discuss 42 Miras that 
have defined the Mira PL-relation and in particular derive pe- 
riods from the MACHO database and compare these to the 
historical period (FGWC, GLE81). In all but three cases they 






Fig. 4. Histograms of the LMC if -band 2MASS magnitude distribution for different "boxes" and /-band amplitude cuts. Shown 
are the histograms for the M-stars (vertical lines), C-stars (hatched), and total. The luminosity function for sequence "D" is likely 
to be incomplete for K < 11 because of the stricter conditions to accept periods above 800 days. 
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Fig. 5. As FigureHfor the SMC. 
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Fig. 7. Colour-Period diagram for the LMC (top left and top right) and SMC (bottom left and bottom right) objects. The period 
with the largest amplitude is plotted. Symbols as in Figure|3] 
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Table 6. First entries in the electronically available table, which list: OGLE-field, OGLE-name, other names, periods (respec- 
tively the historical period-in the majority of cases from the Hughes list or WBP-then the present-day period from OGLE and 
finally other available periods with reference between parentheses), spectral type and references and comments on period changes 



Field 



OGLE-name 



Other Name 



Period SpT Comment 



lmc_scl5 
lmc_scl5 
lmc_scl5 
lmc_scl5 
lmc_scl5 
lmc_scl5 
lmc_scl5 
lmc_scl5 
lmc_scl5 
lmc_scl5 
lmc_scl5 
lmc_scl5 
lmc_scl5 
lmc_scl4 
lmc_scl4 
lmc_scl4 
lmc_scl4 
lmc_scl4 
lmc_scl4 
lmc_scl4 
lmc_scl4 
lmc_scl4 
lmc_scl4 



OGLE050024. 
OGLE050025. 
OGLE050033. 
OGLE050034. 
OGLE050042. 
OGLE050045. 
OGLE050045. 
OGLE050122. 
OGLE05013L 
OGLE050137. 
OGLE050141. 
OGLE050203. 
OGLE050213. 
OGLE050228. 
OGLE050258. 
OGLE050303. 
OGLE050328. 
OGLE050331. 
OGLE050342. 
OGLE050346. 
OGLE050358. 
OGLE050410. 
OGLE050410. 



75-685304.8 
07-690558.0 
74-693017.2 
09-692313.2 
90-690435.2 
01-691422.4 
80-690921.8 
62-684324.0 
98-684130.6 
46-684837.8 

86- 692020.7 
62-692702.9 

87- 692117.6 
98-692010.0 
41-690026.9 
05-691221.7 
77-692747.2 
05-684831.8 
24-691512.4 
90-685753.8 
30-692057.5 
11-692135.2 
21-690743.7 



SHV050038 
SHV050040 
SHV050052 
SHV050052 
SHV050058 
SHV050101 
SHV050101 
SHV050135 
SHV050144 
SHV050150 
SHV050159 
SHV050222 
SHV050231 
SHV050246 
SHV050313. 
KDM.1631, 
SHV050348 
KDM.1683, 
SHV050359 
SHV050401 
SHV050416 
SHV050428 
SHV050426 



5- 685722 
7-691015 
9-693434 

2- 692730 

3- 690851 
9-691838 
9-691338 

1- 684737 

2- 684542 

7- 685250 

8- 692433 

6- 693113 

9- 692527 
9-692418 

7- 690433 

SHV050319.9-691628 
1-693152 

SHV050344.6-685236 

5- 691915 

8- 690157 

6- 692500 
6-692536 
6-691145 



187, 376 
332, 344 
391,432 
315, 310 
162, 163 
317, 333 
276, 283 
231, 243 
401, 386 
146, 146 
386, 360 
255, 246 
180, 177 
310, 312 
307, 302 

1003, 1169 
464, 441 
393, 381 
234, 234 
418, 375 
430, 205 
386, 395 

1353, 1475 
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Fig. 6. Amplitudes versus periods for the LMC (left) and SMC 
(right) variables. Symbols as in Figure|5] 



Table 5. Theoretical colours of M- and C-stars of solar metal- 
licity. Columns 3 and 4 represent the colours for the photo- 
sphere, while, for some models, columns 5 and 6 indicate the 
colours when the star is obscured by circumstellar dust equiva- 
lent to a dimming in K by 0.05 mag. 



SpT 


T e g 


(I -J) 


(J-K) 


(I -J) 


(J-K) 


M0 


3850 


1.18 


1.08 






M3 


3550 


1.39 


1.18 


1.53 


1.34 


M5 


3397 


1.86 


1.27 






M7 


3129 


3.11 


1.37 


3.73 


1.97 


M8 


2890 


3.71 


1.37 






M10 


2500 


4.13 


1.35 


4.45 


1.74 


C 


3600 


1.43 


1.17 






C 


3400 


1.56 


1.30 






c 


3200 


1.72 


1.43 






c 


2800 


2.00 


1.61 






c 


2650 


2.09 


1.61 


2.48 


2.09 



found the periods to be constant over 2 to 3 decades. The three 
are WBP-30 which changed from a 400d Mira into a small am- 
plitude variable with 183d period, GR0537-6740 whose period 
changed from 418 to 367d, and GR17 whose period changed 
from 780 to 729d. In the other objects the change in period was 
less than 3% of the period. 

Figure [8] shows that if-band PL-relation for the stars that 
changed the period by more then 0.1 dex, and with histori- 
cal and recent period connected by a thin line. Several objects 
seem to have changed period considerably and changed pulsa- 
tion mode. The data on these stars have been carefully checked 
in terms of positional association, period determination, am- 
plitude, / and JK magnitudes when available, etc. These stars 
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Fig. 8. if-band Pi-relation for the (LMC) stars that have un- 
dergone a period change over the past two decades of at least 
0. 1 dex. Historical and current period are connected by a line, 
with the current period marked by a solid circle. 

have been marked in Table|6] as well as other stars that changed 
period by more than 10% (30 objects). These smaller changes 
may very well be real but in most cases too small to be re- 
lated to a change in pulsation mode. Most of the changes occur 
from box "C" to "B+" (15 out of 36), and from "C" to "D" 
(8/36). Unfortunately, the original data points of the stars from 
the Hughes list are not available to directly phase the old data 
points with the present day period to verify the change in pe- 
riod. Figure [9] shows some of the lightcurves of the stars that 
possibly changed pulsation mode (complete figure available in 
the electronic edition). As is evident many of them are far from 
being regular. 

Possibly the best studied case of mode switching in a 
Galactic Mira-like variable is that of R Dor. Bedding et al. 
(1998) show that the star switches back and forth between two 
pulsation period of 332 and about 175 days, on a time scale of 
about 1000 days. The star has an accurate (Hipparcos) paral- 
lax and hence can be placed in Figure [8] For any reasonable 
distance to the LMC it implies the star moves back and forth 
between sequences "C" and "B+". 

Figure ^| shows the distribution of the ratio (current pe- 
riod / historical period) near a value of 1 . The observed mean 
is 1.00018. Also a Gaussian fit to this distribution is shown to 
guide the eye, which has a dispersion of 0.03 1 . As stars evolve 
to cooler temperature and lower mass on the AGB one might 
expect an increase of the pulsation period. To verify this the- 
oretically, the synthetic AGB evolution code based on the for- 
mulation in Wagenhuber & Groenewegen (1998) was used for 
a typical star with 2 M Q initial mass and metallicity Z = 0.008. 
The model takes into account the changes in stellar parameters 
over the thermal pulse cycles, and assumes that pulsation takes 
place over the entire TP-AGB, i.e. no cessation of pulsation or 
mode switches during any phase of the stars' TP-AGB evolu- 
tion. For every time step the fundamental and first overtone pul- 
sation period was calculated from P = 0.00851P 194 M- - 90 
and Pi = 0.038P 1 5 M~ - 5 (Wood 1990; R being the stellar 
radius in solar radii, and M the total mass in solar masses). 
Figure [^shows the distribution function of Pq and Pi over a 



timespan of 17 years for such an object, both on a linear and 
logarithmic scale. The AGB evolution is such that on average 
a star evolves to longer period with time, but almost not mea- 
surably so over 17 years. Indeed the averages are 1.00024 for 
fundamental mode and 1.00015 for first overtone. This implies 
that for any individual star one must be able to determine pe- 
riods to fractions of a day to be able to detect period changes 
due to evolution, which seems unrealistic since the lightcurves 
are rather complicated and not mono-periodic in the majority 
of cases. The fact that the observed mean of the period ratios 
over a 17 year timespan is close to the predicted one probably 
indicates good fortune rather than to indicate that, in a statis- 
tical sense, the predicted evolutionary effect has indeed been 
detected. The width of the observed distribution is much wider 
than predicted by the models, and is also wider than expected 
from the errors in the observed period determinations alone. 
The median period of the stars listed in Table [6] is about 290 
days, and in the overall majority of cases such a period has been 
determined with an accuracy of 0.7 days, or better. This implies 
one would expect ratios near unity with an error of about 0.005 
or smaller. This would suggest that the width of the distribution 
is real. 

In addition, there exist a few examples of LPVs whose pe- 
riod decrease over a longer timescale (Wood & Zarro 1981). 
Zijlstra et al. (2002) show that for the Mira R Hya P is about 
-1.6 10~ 3 between 1770 and 1950 AD. Although its period 
has actually stabilised at 385d since then, let us assume it would 
follow its historical trend. Over a 17 year period it would de- 
crease its period by 9.9 days and would have a period ratio, as 
defined and used in Figures [lOlandGj] of 0.975. These values 
are predicted by the stellar evolutionary calculations mentioned 
above (and by the calculations in Wood & Zarro) but with a 
very small probability. The fact that the observed distribution 
in Figure [H)]is wider than expected, based on the errors alone, 
suggests that other phenomena than the "global" increase in pe- 
riod over an AGB stars lifetime dominate this distribution, or 
that a model assumption is incorrect. Either pulsation does not 
take place at all phases of AGB evolution which would influ- 
ence the theoretically predicted distribution function, or other 
physical phenomena play a role. Zijlstra et al. (2002) mention 
and refer to weak chaotic behaviour and the effects of the non- 
linearity of the pulsation in the case of R Hya. 

4.5. Very red objects 

Since the spectroscopically selected sample, for given lumi- 
nosity, will be biased against stars with heavy mass loss, and 
hence fainter / magnitudes, the present section discusses an in- 
frared selected sample with red infrared colours, which will 
complement the sample of spectroscopically selected AGB 
stars. Figure fL2l shows the DENIS Kq — (I — K)q magnitude- 
colour, (/ — J)q — (J — K)q colour-colour, and the 2MASS 
Kq — (J— K)o magnitude-colour, (J— H)q — (H— K)q colour- 
colour diagrams for both LMC and SMC for the spectroscop- 
ically selected sample in the two leftmost columns. Based on 
this diagram it was decided to investigate the pulsation char- 
acteristics of objects that have DENIS (7 - K) > 4.0 and 
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Fig. 9. Representative lightcurves of stars that have changed pulsation mode over the last 2 decades. The fit is indicated by the 
(red) solid line. Crosses indicate data points not included in the fit. Complete figure is available in the electronic form. 



0? + <rk) /{I - K )o < 0.01 or 2mass (J - K) > 2.0 and 
(<t 2 t+v k )/(J-K)o<0M. 

The infrared selected sample contains 577 objects (137 
SMC, 442 LMC) that fullfill these limits and Figure [T3l shows 
the lightcurves of the reddest stars in (J — K) among them 
(the full figure is available in the electronic edition). Tables 
and|8]show the first entries with information similar to Tables|2| 
and [5] The colour-colour diagrams are shown in the two right 
side columns in Figure ^] Some stars from the spectroscopi- 
cally selected sample also appear in the infrared selected sam- 
ple, mostly M-stars for which it is was shown before that the 
later spectral types are reasonably red in (I — K), certainly 
compared to carbon stars for a given (J — K). 

Figure [21 finally shows the if-band PL-relation for the 
large amplitude variables in the infrared selected sample. 



Compared to stars that are spectroscopically known they ap- 
pear to be fainter at a given period. This must be largely due to 
the dust obscuration. An identical effect was shown by Wood 
(2003) who specifically looked at MSX sources. To illustrate 
this, Figure[2]shows the spectral energy distribution of one of 
the stars in Figured namely OGLE050854.2 1-690046.4 or 
MSX 83. For a distance of 50.1 kpc, a luminosity of 22100 L 
is derived and a mass loss rate of 7.6 x 10 -6 M© yr _1 . 

5. Summary and discussion 

This paper addresses several aspects of the pulsational charac- 
ter of late-type stars in the Magellanic Clouds. The main focus 
is on the if-band PL-relation of almost 2300 spectroscopically 
confirmed M-,S- and C-stars. This sample avoids to make the 
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Fig. 11. Distribution function of the ratio of the fundamental (top panels) and first overtone (bottom panels) pulsation period 
over a 17 year timespan, on a linear (left-hand side), and logarithmic (right-hand side) scale for a 2 M Q star with typical LMC 
metallicity. The AGB evolution is such that on average a star evolves to longer periods with time, but not measurably so over 17 
years. 



Table 7. The infrared selected sample. As Table|2] 



Field 

smcjscl 
smcjscl 
smcjscl 
smc_scl 
smcjscl 
smcjscl 
smc_sc2 
smc_sc2 
smc_sc2 
smc_sc2 
smc_sc2 
smcjsc2 
smc_sc2 
smc_sc2 
smc_sc2 
smc_sc2 
smc_sc2 
smc_sc2 
smcjsc3 
smcjsc3 
smc_sc3 
smc_sc3 
smc_sc3 



OGLE- 
OGLE003633 
OGLE003720 
OGLE003732 
OGLE003S01 
OGLE003S30 
OGLE003848 
OGLE004010 
OGLE004014 
OGLE004016 
OGLE004027 
OGLE004032 
OGLE004035 
OGLE004114 
OGLE004116 
OGLE004142 
OGLE004152 
OGLE004156 
OGLE004217 
OGLE004235 
OGLE004240. 
OGLE004243 
OGLE004253 
OGLE004309 



name 

20-735432.5 
78-730944.4 
32-733021.3 
87-734725.7 

92- 733113.6 
77-735601.1 
23-730039.4 
23-724959.2 
25-730115.1 
57-732611.0 

93- 732839.7 
31-731928.7 

43- 731008.8 
95-725216.4 

44- 725853.2 

05- 730827.4 
80-730728.9 
65-733419.1 
80-730920.9 
89-725705.1 

06- 731841.7 
55-725622.2 
55-730922.2 



Period 
0.0 ± 0.000 

128.3 ±0.316 

213.6 ±3.556 
315.0 ±0.508 

100.0 ± 0.026 
280.2 ±0.153 
188.2 ±0.195 

236.1 ± 1.436 

157.5 ± 0.755 

141.4 ± 0.126 

551.2 ±0.643 

508.7 ± 0.765 

625.2 ± 3.056 

384.3 ± 0.446 

415.8 ± 0.343 
366.7 ±0.185 

241.5 ±0.218 

131.9 ±0.026 

121.6 ±0.159 

221.6 ± 0.295 
382.1 ±0.126 
359.3 ± 1.282 
383.6 ±0.215 



Amp 
0.000 
0.117 
0.025 
0.142 
0.243 
0.588 
0.128 
0.206 
0.054 
0.086 
1.263 
0.881 
0.943 
0.629 
0.336 
0.580 
0.125 
0.489 
0.066 
0.173 
0.790 
0.539 
0.873 



Period 
0.0 ± 0.000 
424.9 ± 0.680 

1192.8 ± 99.000 
494.8 ±0.189 
158.7 ±0.052 
640.1 ±9.035 
351.5 ±0.151 

250.3 ± 1.534 
480.5 ± 0.925 
269.5 ±0.135 

0.0 ± 0.000 
0.0 ± 0.000 
862.6 ± 30.639 
1484.6 ± 16.917 
0.0 ± 0.000 
1538.7 ± 6.836 

475.4 ± 0.248 
0.0 ± 0.000 

493.7 ± 0.244 

491.8 ±0.870 
1 127.0 ± 4. 1 IS 

1845.9 ± 99.000 
0.0 ± 0.000 



Amp 
0.000 
0.734 
0.597 
0.640 
0.380 
0.179 
0.774 
0.219 
0.853 
0.481 
0.000 
0.000 
0.352 
0.688 
0.000 
0.726 
0.331 
0.000 
0.899 
0.681 
0.245 
0.431 
0.000 



Period 
0.0 ± 0.000 
0.0 ± 0.000 
2056.4 ± 99.000 
0.0 ± 0.000 
384.2 ± 1.362 
1499.2 ± 16.100 
1979.9 ± 24.789 
1095.9 ± 3.42(1 
1503.1 ± 19.637 
515.5 ± 1.259 
0.0 ± 0.000 
0.0 ± 0.000 
0.0 ± 0.000 

548.0 ± 8.036 
0.0 ± 0.000 
0.0 ± 0.000 

935.1 ± 3.473 
0.0 ± 0.000 

1959.7 ± 27.705 
1395.1 ±6.272 
0.0 ± 0.000 
0.0 ± 0.000 
0.0 ± 0.000 



Amp 
0.000 
0.000 
0.736 
0.000 
0.079 
1.448 
0.510 
0.203 
0.648 
0.176 
0.000 
0.000 
0.000 
0.121 
0.000 
0.000 
0.168 
0.000 
0.416 
0.861 
0.000 
0.000 
0.000 



16.61 
19.04 
18.44 
15.58 
15.31 
17.22 
15.75 
14.51 
18.10 
15.76 
18.66 
18.85 
18.74 
15.92 
15.85 
1 6.69 
14.21 
15.61 
17.44 
17.24 
17.04 
15.89 
I S.43 



16.11 ±0.03 
99.00 ± 9.99 
99.00 ± 9.99 
16.19 ±0.05 
14.74 ± 0.02 
16.76 ± 0.05 
15.78 ±0.05 
14.80 ± 0.02 
17.00 ±0.11 
15.86 ± 0.04 
17.30 ±0.19 
99.00 ± 9.99 
99.00 ± 9.99 
16.18 ±0.05 
16.42 ±0.07 
99.00 ± 9.99 
14.09 ±0.01 
15.04 ±0.01 
17.49 ±0.18 
16.53 ± 0.08 
99.00 ± 9.99 
15.76 ± 0.04 
99.00 ± 9.99 



Jdeni, 

13.90 ± 0.03 
99.00 ± 9.99 
99.00 ± 9.99 
13.79 ± 0.03 
13.00 ± 0.01 
15.18 ±0.17 
13.68 ± 0.04 
12.66 ± 0.02 
15.87 ± 0.25 
14.08 ± 0.04 
14.20 ± 0.07 
99.00 ± 9.99 
99.00 ± 9.99 
13.85 ± 0.04 
13.52 ± 0.03 
99.00 ± 9.99 
12.13 ± 0.01 
13.24 ± 0.01 
14.66 ± 0.09 
15.24 ±0.11 
99.00 ± 9.99 
13.61 ± 0.03 
99.00 ± 9.99 



A'deni, 

11.00 ±0.02 
99.00 ± 9.99 
99.00 ± 9.99 
99.00 ± 9.99 
10.54 ± 0.01 
12.17 ±0.04 
10.90 ± 0.03 
10.52 ± 0.02 
12.26 ±0.10 

11.74 ±0.06 
10.79 ± 0.03 
99.00 ± 9.99 
99.00 ± 9.99 
11.19 ±0.03 

11.75 ±0.05 
99.00 ± 9.99 
10.10 ±0.01 
10.78 ± 0.01 
11.28 ±0.04 
11.40 ±0.03 
99.00 ± 9.99 
10.97 ± 0.02 
99.00 ± 9.99 



JgMASS 
14.14 ±0.03 
16.04 ±0.11 
16.08 ± 0.09 
13.81 ±0.02 
13.08 ± 0.02 
15.20 ± 0.07 
14.43 ± 0.03 
12.58 ± 0.03 
14.28 ± 0.04 
13.38 ±0.03 

15.61 ± 0.08 
16.28 ±0.12 
16.26 ±0.11 
14.00 ± 0.03 
13.67 ±0.02 
13.90 ±0.02 
12.28 ± 0.02 
13.60 ± 0.03 
15.30 ± 0.05 
13.22 ±0.03 

13.62 ±0.05 
13.13 ±0.02 
14.18 ±0.07 



#2MASS 

12.55 ± 0.03 
13.63 ± 0.03 
14.62 ± 0.05 
12.23 ± 0.02 
11.68 ±0.03 

13.56 ± 0.05 
12.80 ± 0.03 
11.30 ±0.02 
12.46 ± 0.03 
12.23 ± 0.02 
13.39 ± 0.04 

14.18 ±0.05 
13.97 ± 0.05 
12.33 ± 0.03 

12.00 ± 0.03 

12.27 ± 0.03 

11.01 ±0.02 

12.19 ± 0.03 
13.29 ± 0.02 
11.61 ±0.03 
12.23 ± 0.05 
11.77 ± 0.02 

12.28 ± 0.05 



^2MASS 

1 1 .25 ± 0.03 
11.85 ±0.02 
13.60 ± 0.05 
10.79 ± 0.02 

10.72 ± 0.02 
12.28 ± 0.03 
11.45 ±0.02 
10.51 ±0.02 
10.97 ± 0.02 
11.28 ±0.02 
11.57 ±0.03 
12.30 ±0.03 

11.97 ±0.03 

10.98 ± 0.02 
10.64 ± 0.02 
11.05 ±0.05 
10.16 ±0.02 
11.02 ± 0.02 

11.73 ±0.03 
10.16 ± 0.02 
11.08 ±0.03 

10.78 ± 0.03 

10.79 ± 0.04 



clearly incorrect approximation made in other studies that M- 
stars and carbon stars can be separated at a colour (J — K) = 
1.4. 

The present observations however do not allow the pre- 
sentation a comprehensive picture of the evolution of pulsa- 
tion periods of M- and C-stars. This would require a more de- 
tailed (AGB star) population synthesis study including pulsa- 
tion properties. As previous MCs studies have found for AGB 
stars in general (Wood et al. (1999), Wood (2000), Noda et 



al. (2002), Lebzelter et al. (2002), Cioni et al. (2003), Ita et 
al. (2004a,b) and Kiss & Bedding (2003, 2004)), it is found 
specifically that both M-, and C-stars tend to occupy preferen- 
tially sequences B+ and then C for increasing amplitude. For 
a given amplitude, the C-stars tend to have the longer period, 
and, for a every sequence they are more luminous. This effect 
has previously been observed in MC clusters (e.g. FMB) and 
is in qualitative agreement with evolutionary calculations that 
predict that C-stars evolve from M-stars. 
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Fig. 12. LMC (top 2 rows) and SMC (bottom 2 rows) colour-colour diagrams using 2MASS and DENIS photometry for the 
spectroscopically selected sample (left 2 columns) and the infrared selected sample (right 2 columns). Note the difference in 
scale! Symbols as in Figure|5] with dots indicating objects without spectroscopic classification. 



Many objects have one period that falls in box "D". In the 
spectroscopically selected sample, 211 of 859 SMC stars (= 
24.6%) have a period that falls in box "D", and 318 of the 1418 
LMC stars (= 22.4%, namely 229/1064 = 21.5% of C-stars and 
89/354 = 25.1% of M-stars). As will be discussed below, for 
at least some stars of the IR selected sample this is due to the 
fact that they are weaker in K because of dust obscuration. 
For the overall majority of the stars in the spectroscopically 
selected sample this is not an issue. The reason why some late- 
type stars appear on that location of the Pi-diagram is unex- 
plained, see the discussion in Olivier & Wood (2003), Wood 
(2003, 2004). The classical large-amplitude Mira variables ap- 
pear on sequence "C" (see the last panels in Figure [3} and are 
believed to be fundamental-mode pulsators, hence longer (ra- 
dial mode pulsation) periods should not exist. The present pa- 
per does not shed light into the nature of the LSP phenomenon, 
except that Figure|4]indicates that the A" -band luminosity func- 



tion for sequence "D" is essentially the same for the M- and 
C-stars, while for sequences "A,B,C" the C-stars are brighter, 
as expected from an evolutionary point of view. The luminos- 
ity function of objects on sequence "D" is fainter from those of 
the other sequences. This is due to the fact that periods longer 
than 800 days are underrepresented because of stricter selection 
rules. This affects predominatly sequence "D" objects with Ko 
<,10.5 mag. The fact that the fraction of all object with peri- 
ods on sequence "D" is essentially the same for C- and M-stars 
and that the if-band luminosity of C- and M-stars of sequence 
"D" objects is very similar suggests that the LSP phenomenon 
is unrelated to chemical type, and hence seems unrelated to a 
pulsation phenomenon. 

For a few hundred variables it was possible to look for pe- 
riod changes over a timespan of typically 17 years. Almost all 
come from the studies by Hughes (1989) and Hughes & Wood 
(1990). They identified medium to large amplitude variables 
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figl3.jpg 



Fig. 13. Lightcurves of the reddest stars in (J — K). Note the faintness in I. Nine of the 15 stars have been detected in the MSX 
survey. Complete figure is available in electronic form. 



from photographic material using typically 21 observations in 
the time span 1977 to 1984. Out of 370 objects, 36 have been 
identified that seemingly changed pulsation mode (or at least 
changed "box") between ~4980 and the time of the OGLE ob- 
servations, and another 30 objects that changed pulsation pe- 
riod by more than 10%. This ratio of about 10% (36/370) is 
similar to the study of GLE03 who found large period changes 
in 3 out of 42 Mira variables studied. A caveat is that the orig- 
inal historical data points have never been published, and it 
would certainly be preferable to be able to phase the old data 
with the current period to see if in fact the period change is 
real. For the moment I consider the 10% change of pulsation 
sequence over ~17 years as an upper limit. By comparison, 
Zijlstra & Bedding (2003) find that only of order of 1% of 
well-known Miras show evidence for period changes. The un- 
derstanding for MC objects may improve because of ongoing 
(e.g. OGLE-III) and future surveys. 



Finally a sample of stars was studied selected on infrared 
colours, namely redder than the majority of the spectroscop- 
ically selected sample. It should be pointed out that there is 
no proof that these are AGB stars (except for the few ones in 
overlap with the spectroscopically selected sample) and they 
make suitable targets for spectroscopic follow-up to determine 
their spectral type. Many of these stars also have a period lo- 
cated in box "D" but in this case the effect of obscuration by 
dust must be considered. This also has implications when using 
samples of variables to determine if -band PL-relations. In the 
last section one explicit example was shown, namely MSX 83, 
for which the SED was constructed and fitted with a dust ra- 
diative transfer model. Its period of 61 1 day, and K magnitude 
of 10.58 would put it in "box D". However, when running the 
radiative transfer model without mass-loss the if -magnitude 
brightens to 9.17 mag (this value is somewhat dependent on 
the central star model atmosphere assumed) putting it on the 
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Fig. 14. AT-band Pi-relation for the large amplitude variables in the IR selected sample, for the SMC (left) and LMC (right). 
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Fig. 10. Distribution function of the ratio of the historical to 
the current period near a value of unity. The mean of all values 
is 1.00179. The solid line is a Gaussian fit with a a of 0.031, 
centered at 0.9981. 

Table 8. The infrared selected sample. As Table[3] 



Field 



OGLE-name 



SpT Comments 



smcscl OGLE003633. 20-735432.5 

smc_scl OGLE003720.78-730944.4 

smcscl OGLE003732.32-733021.3 

smcscl OGLE003801.87-734725.7 

smcscl OGLE003830.92-733113.6 

smcscl OGLE003848.77-735601.1 

smcsc2 OGLE004010.23-730039.4 

smcsc2 OGLE004014.23-724959.2 

smcsc2 OGLE004016.25-730115.1 

smcsc2 OGLE004027.57-732611.0 

smcsc2 OGLE004032.93-732839.7 

smcsc2 OGLE004035. 31-731928.7 

smcsc2 OGLE004114.43-731008.8 

smcsc2 OGLE0041 16.95-725216.4 

smcsc2 OGLE004142.44-725853.2 

smcsc2 OGLE004152.05-730827.4 

smcsc2 OGLE004156.80-730728.9 

smcsc2 OGLE004217.65-733419.1 

smcsc3 OGLE004235. 80-730920.9 

smcsc3 OGLE004240.89-725705.1 

smcsc3 OGLE004243.06-731841.7 

smcsc3 OGLE004253. 55-725622.2 

smcsc3 OGLE004309.55-730922.2 



MH.0400 



[WSH2003] a- 8 




log A [^m] 

Fig. 15. SED of one of the stars in the IR selected sample, MSX 
83 annex OGLE050854.2 1-690046.4. Plotted are the ogle /, 
2MASS JHK, MSX A band and IRAS 12 and 25 /jm fluxes. 
A carbon star central star model is assumed. 

extension of box "C", consistent with the expected location for 
a star with an pulsation amplitude of 0.66 mag. 

In reverse this implies that when studying the A"-band PL- 
relation, and when multi-colour data is available, a cut-off in 
colour should be applied in order to avoid a bias by including 
stars that are dust obscured in K. Although this seems obvi- 
ous, a quantification of where this cut-off should be placed and 
its actual application are rare in the literature; In fact I could 
only find one instance. Glass et al. (1995) mention they exclude 
some faint outliers with (K — i) ~ 2 in the Sgr I bulge field 
(which corresponds to roughly (J — K) ~ 3.5, Glass, 1986), 
but did not impose a colour criterium a-priori. In other papers 
where dust obscuration in the studied variables should play a 
role the bolometric Pi-relation is studied (e.g. WFLZ), which 
circumvents this problem in a natural way. 

It is difficult to give an exact colour cut-off to apply, since 
this depends on the colours involved, the dust properties and the 
evolution of mass-loss on the AGB. Based on the calculations 
presented in the last columns in Table|5]stars with colours ( J — 
K) > 2.0 should certainly be avoided, and a stricter criterium 
would be to include M-stars only when (J — K)o < 1.4 and 
C-stars only when (J — K)q < 1.7. 

Finally, applying these latter cut-offs to the spectroscopi- 
cally selected stars in box "C" with amplitudes > 0.45 mag 
and > 0.15 mag, to improve the statistics, the AT -band PL- 
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relations listed in Table|9]have been derived, where also some 
relations from the literature are listed. The period distribution 
of these samples is shown in Figure^] It would be interesting 
to fit the SEDs of all infrared selected stars to be able to include 
the dust-correced A'-magnitudes in these PL-relations. 

There is very good agreement between the new relations 
(including what appears to be the first Mira LT-band PL- 
relation in the SMC) and previous works. The formal error on 
the zero point and slope have become smaller, because of the 
larger sample size, but the overall rms are still larger because 
the photometry used in FGWC and GLE03 are averages over 
the lightcurve while the present data are at best averages of two 
measurements. 
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Fig. 16. Period distribution of the LMC (left panel) and SMC (right panel) variables in box "C" with /-band amplitudes larger 
than 0. 15 magnitudes. Shown are the histograms for the M-stars (vertical lines), C-stars (hatched), and total. 
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Appendix A: Step 1 : Fourier analysis and light 
curve fitting 

At the heart of the first code are the NUMERICAL RECIPES subrou- 
tines fasper and mrqmin that perform a Fourier transformation and a 
(weighted) linear least-squares fitting, respectively. The function that 
is fitted to the data has the following form: 

n max 

I(t) = I + Y ( A i sin(27r t e h ) + Bi cos(2tt t e fi )j (A.l) 

This is a more suitable form for fitting purposes than the equivalent 
form (with cu = e f ): 

'-Umax 

I(t) = J + ^2 Ci sin(27r t Wi + (pi) 
»=i 

Up to three periods are fitted (n max = 3), following Wood et al. (1999). 
These two subroutines are called alternatingly, subtracting the best fit 
sofar from the data, and then performing a Fourier analysis on the 
residual. If a significant peak is found, a fit including the next term 
in the series in Eq. 1 is included. This is repeated until no significant 
peak in the power spectrum is found or the third period has been fitted. 

The Fourier analysis is done with the subroutine fasper. Inputs to 
it are the time and magnitude arrays. In addition one has to specify 
two parameters, ofac and hifac, that indicate a "typical oversampling 
factor" and the maximum frequency in terms of a "typical" Nyquist 
frequency. 

In the present work ofac = 22 and hifac = 0.8 are used. The latter 
parameter is the determining factor in both the computational speed, 
and the shortest period that can be found. For example, some tests 
to correctly identify the main period in an RR Lyrae object with a 
known period of 0.55 days required hifac = 10.0. In this configuration 
the code would be more than a factor of 10 slower. At the same time 
it implies that in the configuration used in the present paper there is 
a bias in the detection of periods shorter than about 6 days, of no 
consequence as the focus is on AGB stars. 

The output of fasper are the frequency where the peak occurs 
and a number indicating a significance. One of the main parameters 
in the code is to provide the critical cut-off above which a period is 
not considered to be significant. In the present work significance = 
5.5 x 10 -11 is used, and this was determined empirically, by visually 
inspecting many lightcurves. 

The code can be run in a single-star mode (for fine tuning) or in an 
automatic mode. It should be pointed out that some of the features and 
parameters just described and that will be described below have been 
determined only empirically. In fact, the cycle of Steps 1 and 2, and 
partly Step 3, has been repeated a few times to check the various steps 
in detail and come to the final choice of the parameters. The content 
of the code is now described in detail. 

- Read in the files with 2MASS and DENIS objects that are within 
3" of OGLE objects. 

Read in the file which contains the path names to the 68 000 files 
which contain the /-band data. 

The steps below are either done for one star (in single-star mode) 
or for all. 



- Read in the path-name (e.g. smc_sc5/OGLE005013.32- 
73 11 12.6. mag), and construct the right ascension and declination 
from it. Read in the individual Julian Dates, /-band magnitudes 
and errors. If there are less than 11 data points, exit. Otherwise 
enter the light curve fitting subroutine. 

- Loop over the data points, and remove those with error "—99". 
Loop over the remaining data points and determine the largest 
time gap between two subsequent data points in time. If there are 
less than 30 data points or a time gap larger than 250 days, exit. 

- Order the data points according to the error bar in the I- 
magnitude. Automatically reject 5% of the data with the largest 
errors. Order the data points according to the /-magnitude. 
Automatically reject the 3 brightest data points. This last step was 
introduced when it became clear that a few very bright (spuri- 
ous ?) outliers in faint sources severely biased the light curve fit- 
ting because these outlying points carry so much weight. In bright 
sources this procedure has no bearing on the outcome. The final 
data arrays with time, /-magnitude and error in /-magnitude are 
constructed, and written to file, so that they can be read in by an in- 
dependent time-analysis code (PERIOD98, Sperl 1998) in single- 
star mode. 

- I. The best fit so far is subtracted from the data. At the first in- 
stance a constant is subtracted. The residual is inputted to fasper. 
If the significance of the peak is above the cut-off, or already n max 
periods have been determined, goto II. 

In single-star mode the frequency sometimes has to be set to twice 
or half the frequency found by fasper at this stage in order to let 
the program converge to the correct result. 

- First guesses for the amplitudes (see Eq. 1) are determined, and 
then the data plus parameters are inputed to mrqmin. At least two 
calls to mrqmin are made, and these calls are then repeated until 
subsequent \ 2 values agree to within 1 % (or until a maximum of 
10 iterations is reached). The \ 2 value is stored. 

- For the improved fit, for each data point the deviation in terms of 
sigma is determined and if the data point with the largest deviation 
deviates more than 10a from the fit, the error bar of that point is 
set to a very large value, so that that point is ignored in the further 
analysis. 

- goto I. 

- II. Plot the light curves including the best fit, and write all relevant 
output to file. 

- Do the final correlation with DENIS and 2MASS using a search 
radius of 1.5". This takes into account the typical accuracy of 
the coordinates in the OGLE, DENIS and 2MASS surveys. Multiple 
matches are allowed for. The relevant data is written to file. 

Appendix B: Step 2: Selecting AGB stars and 
general statistics 

In the second code a preliminary list of LPV candidates is determined, 
by applying (if one wishes) selections in magnitude, periods and am- 
plitudes and by eliminating known non-LPVs and correlating with 
known LPVs and/or AGB stars. These issues are discussed here. 

B. 1. Known sources 

The advantage of compiling a list of variable objects in the direction 
of the MCs that are known not to be LPVs is twofold. First, it imme- 
diately limits the number of sources to be inspected visually in Step 3. 
Even more importantly, knowing the light curves of known non-LPVs 
helps identifying other such kind of objects in Step 3. 

The literature was scanned for lists of objects identified in OGLE 
and other microlensing surveys MACHO, MOA, EROS in the direction 
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of the MCs (hence not restricted to the OGLE fields). Eclipsing bi- 
naries (EBs) were included from Udalski et al. (1998; OGLE, 1459 
in SMC), Wyrzykowski et al. (2003; OGLE, 2580 in LMC), Alcock 
et al. (1997; MACHO, 637 in LMC), Bayne et al. (2002; MOA, 167 
in SMC), Grison et al. (1995; EROS, 79 in LMC). RV Tau objects 
from Alcock et al. (1998; MACHO, 33 in LMC). R CrB stars from 
Alcock et al. (2001; MACHO, 17 in LMC). A few known and many 
new candidate QSO from Eyer (2002; OGLE, 133 in SMC plus LMC), 
and Geha et al. (2003; MACHO, 59 SMC + LMC). Blue variable ob- 
jects (possibly related to the Be phenomenon) from Mennickent et al. 
(2002; OGLE, 1056 in SMC), Mennickent et al. (2003; OGLE, 30 in 
SMC plus LMC), Eyer (2002; OGLE, 36 in SMC plus LMC), Keller 
et al. (2002; MACHO, 1280 in LMC). RR Lyrae stars from Soszynski 
et al. (2002, OGLE, 556 in SMC; 2003, 7612 in LMC) and Alcock 
et al. (2000; MACHO, 283 in LMC). Finally, Cepheids from Afonso 
et al. (2003; EROS, 880 in SMC and LMC), second overtone (SO) 
cepheids from Alcock et al. (1999; MACHO, 47 in LMC) and funda- 
mental mode (FU), first overtone (FO), SO and double-mode (FU/FO 
and FO/SO) cepheids from Udalksi et al. (1999a,b,c,d; OGLE, 3492 
in SMC + LMC). 

This amounts to a total of 20863 objects (including double entries 
coming from different sources). 

In Step 3 the SIMBAD database is queried and so previous obser- 
vations will normally be identified in this way. On the other hand, not 
all (recent) data is yet included there. Therefore a list of known or sus- 
pected LPVs, mostly from recent large survey work, was compiled. 

Two of the largest 'old' surveys for LPVs are those by Hughes 
(1989) and Reid, Hughes and Glass (1995). I was not able to trace the 
relevant tables in these papers in computer readable form, but through 
P. Wood I obtained a table originally prepared by S. Hughes that seems 
to combine the data from these two papers. Eliminating double entries 
it contains information on 1317 LPVs in the LMC (I will refer to this 
list as the Hughes-list from now on). 

From the MOA survey the 146 LPVs in the LMC discussed by 
Noda et al. (2002) are considered. 

M. Schultheis kindly made available the unpublished periods, am- 
plitudes, magnitudes of the about 470 LPVs in the LMC bar discov- 
ered by the AGAPEROS survey and discussed in detail by Lebzelter et 
al. (2002) 

From Cioni et al. (2003) the 458 objects in the SMC are consid- 
ered that have been detected in a mini-survey with ISOCAM (Loup et 
al. 2004, in prep.), have reliable DENIS and 2MASS counterparts (in 
the 2nd incremental data release) and have a MACHO lightcurve. 

P. Wood kindly made available the MACHO id numbers, coordi- 
nates, periods and photometry of the 1560 objects he studied in Wood 
et al. (1999) and Wood (2000) in the bar of the LMC. 

From Feast et al. (1989), using additional information provided in 
Glass & Lloyd-Evans (2003) and references therein, 58 LMC objects 
are considered that de facto have been defined the Mira PL-relation. 

Finally, a list of 45 IRAS sources was compiled from Wood et al. 
(1992), Wood (1998), Whitelock et al. (2003), that have been moni- 
tored in the infrared and were found to have well determined periods. 

The total list of known LPVs contains 4053 objects. 

In addition a list of spectroscopically confirmed M-, S-, and C- 
stars and supergiants with accurate coordinates (listed to 1" or better) 
was compiled from lists in Westerlund et al. (1981; WOH), Prevot 
et al. (1983; PMMR), Wood et al. (1985; WBP), Sanduleak (1989; 
SkKM), Reiberot et al. (1993; RAW), Kontizas et al. (2001; KDM), 
Morgan & Hatzidimitriou (1995; MH), Demers et al. (1993), Kunkel 
et al. (1997), Groenewegen & Blommaert (1998; GB98), Kunkel et al. 
(2000), Loup et al. (2003, specifically dealing with the Blanco fields 
as discussed in Blanco et al. (1980; BMB) and Frogel & Blanco 1990) 



and Cioni et al. (2001), for a total of 13175 stars (including some 
double-entries). 

As for some of the stars in the Hughes-list a spectral type has been 
determined, in fact 13451 stars (including double entries) in this list 
have a spectral type assigned. Correlating this list with itself, using a 
4" search radius, reveals 12631 unique entries (2899 C, 19 M, S in 
the SMC, 8117 C, 1580 M, 16 S in the LMC). This list is strongly 
biased towards carbon stars because of the very large surveys in both 
LMC and SMC (notably RAW and KDM), and the lack of similar 
surveys for M-stars. 

B.2. Cuts in magnitude, amplitude and period 

In the present paper no cuts on magnitude, amplitude or period have 
been made, as the sample discussed below will be restricted to spectro- 
scopically confirmed M-, S-, C-stars. For reference, in Groenewegen 
(2004) the following cuts were applied to the SMC data to select the 
LPV candidates: mean OGLE I < 16.8 mag, any of the fitted periods 
> 50 days, and any of the fitted amplitudes (in the classical sense, 
C t = sJAf + Bf from Eq. 1) > 0.05 mag. 

In hindsight, it would have been preferable to have had a faint 
limit (maybe at / ~ 18-18.5) imposed as quite a few spectroscopi- 
cally confirmed M,S,C-stars were positionally matched with very faint 
(mean / £19) OGLE objects which were often barely variable, and 
clearly not the counterpart of the late-type stars. 

B.3. More details 

The content of the code is now described in detail. 

- Read in the output file of Step 1 , and the files with the coordinates 
of the known non-LPVs, and the known LPVs and spectroscopi- 
cally confirmed M-,S-,C-stars. 

- Only periods with an accuracy better than 2% are retained. 

- The selection on mean /-magnitude, amplitude and period may be 
carried out. In the present paper no cuts have been applied. 

- Objects that appear twice in different OGLE fields are identified 
and the one with the lowest x as determined in Step 1 is retained. 

- The remaining objects are correlated on position (search radius 
2") with the known non-LPVs. The identifiers of both the known 
non-LPVs and the remaining objects are written to file. This sep- 
aration assumes that the classification in the literature of an object 
as Cepheid, RR Lyrae, eclipsing binary, etc, is correct. 

- The remaining objects are correlated on position (search radius 
4") with the known LPVs and spectroscopically confirmed M-,S- 
,C-stars. The identifiers of the cross-correlated objects are written 
to file in Latex format. 

In fact, to allow for small differences in the astrometry, corrections 
have been determined and applied, as discussed later. 

- Independent of the selection of potential LPVs, in the second part 
of the code, general properties are investigated, like the positional 
match between OGLE and 2MASS, and OGLE and DENIS sources, 
and the difference between 7 og ie and /denis, </2mass and Jdenis, 
and K2ma.ee and K dcnis . 

Appendix C: Step 3: Visual inspection and 
literature study 

This last step is time consuming and involves several checks. 

- A visual inspection of the fit to the data. If needed the fit is re-done 
with the code in single-star mode, and usually involves "twisting" 
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the first frequency it finds to half or twice its value. On some oc- 
casions the data is analysed with an independent code, PERIOD98 
(Sperl 1998). 

- Cross-correlation with the SIMBAD database. The coordinate list 
with candidate LPVs is send to the batch queue server of SIMBAD. 

- Checking and completing of automatically generated tabular ma- 
terial in Latex format. 
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